Corals adapted to extreme and fluctuating seawater pH increase calcification rates and have unique symbiont communities

Abstract Ocean acidification (OA) is a severe threat to coral reefs mainly by reducing their calcification rate. Identifying the resilience factors of corals to decreasing seawater pH is of paramount importance to predict the survivability of coral reefs in the future. This study compared corals adapted to variable pHT (i.e., 7.23–8.06) from the semi‐enclosed lagoon of Bouraké, New Caledonia, to corals adapted to more stable seawater pHT (i.e., 7.90–8.18). In a 100‐day aquarium experiment, we examined the physiological response and genetic diversity of Symbiodiniaceae from three coral species (Acropora tenuis, Montipora digitata, and Porites sp.) from both sites under three stable pHNBS conditions (8.11, 7.76, 7.54) and one fluctuating pHNBS regime (between 7.56 and 8.07). Bouraké corals consistently exhibited higher growth rates than corals from the stable pH environment. Interestingly, A. tenuis from Bouraké showed the highest growth rate under the 7.76 pHNBS condition, whereas for M. digitata, and Porites sp. from Bouraké, growth was highest under the fluctuating regime and the 8.11 pHNBS conditions, respectively. While OA generally decreased coral calcification by ca. 16%, Bouraké corals showed higher growth rates than corals from the stable pH environment (21% increase for A. tenuis to 93% for M. digitata, with all pH conditions pooled). This superior performance coincided with divergent symbiont communities that were more homogenous for Bouraké corals. Corals adapted to variable pH conditions appear to have a better capacity to calcify under reduced pH compared to corals native to more stable pH condition. This response was not gained by corals from the more stable environment exposed to variable pH during the 100‐day experiment, suggesting that long‐term exposure to pH fluctuations and/or differences in symbiont communities benefit calcification under OA.


| INTRODUC TI ON
Coral reefs have been suggested to be at risk of disappearance in the coming decades due to their sensitivity to the deadly trio of stressors including ocean warming (OW), acidification (OA), and deoxygenation (OD) (e.g., Hoegh-Guldberg et al., 2017;Hughes et al., 2018Hughes et al., , 2020. These stressors are expected to increase over the coming decades, leading to a dramatic decline of coral reefs by as early as 2050 (Hoegh-Guldberg et al., 2017;IPCC Report, 2021).
More generally, coral resilience to extreme conditions includes trade-offs in the main physiological traits, such as reduced calcification rate (Bay & Palumbi, 2016;Camp et al., 2020;Cunning et al., 2015) and photosynthesis (Camp et al., 2020), increased photodamage (Silverstein et al., 2017), and change in the autotrophic nutrition capacities (Wall et al., 2020). For instance, corals with thermotolerant symbionts have been reported to resist bleaching but grow more slowly (Jones & Berkelmans, 2010). Most of these trade-offs remain poorly tested on corals already adapted to conditions close to or even worse than the ones predicted for the end of this century. Yet, extreme environments, such as mangrove lagoons have been identified to be faithful natural laboratories for studying species living under the combination of future extreme conditions in OA, OD, and OW .
The semi-enclosed lagoon of Bouraké in New Caledonia is one of these extreme environments which has been recently environmentally and ecologically characterized (Camp et al., 2017;Maggioni et al., 2021;Tanvet et al., 2022). Bouraké has fluctuating acidified, warm, and deoxygenated conditions, although a highly-diversified coral reef develops (Camp et al., 2017). The environmental variability in the dissolved oxygen, temperature, and pH in Bouraké is directly related to the tidal cycle with changes on a single day by up to 4.91 mg O 2 L −1 , 6.50°C, and 0.69 pH T units (Maggioni et al., 2021), respectively. Such environmental conditions likely enhanced coral respiration rates and total chlorophyll content, while reducing photosynthesis and calcification rates (Camp et al., 2017;Jacquemont et al., 2022). In addition to the particular physiological response of some corals from Bouraké, Camp et al. (2020) also found differences in host-microorganisms associations (Symbiodiniaceae and bacteria) between Bouraké and an adjacent reef. Although species-specific, the differences in the Symbiodiniaceae found in three coral species (Camp et al., 2020) highlight the potential role that Symbiodiniaceae play in influencing coral holobiont resilience to environmental stress (Sampayo et al., 2008). Corals in Bouraké have likely either adapted to, or sustain key trade-offs to achieve resilience in an environment somewhat reproducing future climate conditions.
In the present study, we conducted a 100-day OA experiment in aquaria using three coral species (Acropora tenuis, Montipora digitata, and Porites sp.) from the highly variable Bouraké and a less variable adjacent reference reef to assess: (i) whether corals adapted to ambient or fluctuating pH conditions alter their metabolic rates and calcification under different levels of OA (i.e.,7.54,7.76,8.11,, and (ii) whether Symbiodiniaceae communities are distinct between habitats and treatments at the end of the experiment.
In addition, we considered both static pH NBS (7.54 ± 0.08; 7.76 ± 0.07; 8.11 ± 0.05) and variable (ranging from 7.56 ± 0.07 to 8.07 ± 0.07) pH NBS conditions to assess the role that fluctuating pH may play in the success of Bouraké corals, and whether corals adapted to stable pH can acclimate when exposed during 100 days to variable pH. We hypothesize that corals from Bouraké exhibit enhanced physiological traits and a distinct Symbiodiniaceae community, compared with corals from the less variable reef, and that the natural diel fluctuations in seawater pH promote their resilience to OA.

T A X O N O M Y C L A S S I F I C A T I O N
Ecophysiology, Evolutionary ecology, Genomics, Global change ecology 2 | MATERIAL S AND ME THODS
Environmental conditions at both sites have been previously characterized across several days in February, March, and May 2016 (Camp et al., 2017) and across numerous surveys (from 22 to 72 semidiurnal tidal cycles, depending on environmental parameters and site) from 2017 to 2020 (Maggioni et al., 2021). In Bouraké, dissolved oxygen (DO) and pH regularly fluctuated according to the tide (ca. 1 m tidal range), from 1.87 to 7.24 mg O 2 L −1 (mean of 5.23 ± 0.89 mg O 2 L −1 recorded during 42 semidiurnal tidal cycles from 2016 to 2020) and from 7.23 to 8.06 pH T (mean of 7.67 ± 0.23 pH T recorded during 72 semidiurnal tidal cycles from 2016 to 2020), respectively (Maggioni et al., 2021). Seawater temperature varied from 17.50°C in winter to 33.80°C during summer, when the temperature might change by up to 6.50°C in a single day (mean of 25.63°C ± 2.85°C during a continuous record from January 2019 to April 2020 with 10-minute interval) (Maggioni et al., 2021). In comparison to Bouraké, DO, pH, and temperature were more stable at the reference site (hereafter reference) during equivalent periods with means of 6.45 ± 0.95 mg O 2 L −1 , 8.02 ± 0.04 pH T , and 25.25°C ± 1.89°C, respectively (Maggioni et al., 2021). For example, at the end of the austral summer from February to March 2016, the reefs in the Bouraké lagoon spent 44% of the time at pH T of 7.7-7.8, and 71% of the time at temperatures predicted for the end of the century under scenario RCP7.0 (Camp et al., 2017). Notwithstanding such extreme and chronic conditions, 66 coral species form reefs all along the Bouraké lagoon (Maggioni et al., 2021).
Fifteen mature colonies (30-60 cm diameter) of the branching A. tenuis (Dana, 1846) andM. digitata (Dana, 1846), and the massive Porites sp. were sampled at each site at ca. 2-meter depth. Only colonies at a distance of at least 5 m from each other were collected to limit the risks of clonality in our sampling. From each branching colony, four terminal portions of branches (ca. 3-5 cm long) were collected using pliers. For Porites sp., four samples (mean of ca. 8 cm 2 ) were collected from each colony using a 3-cm-diameter steel tube and a hammer. Coral fragments were transported in a cooler to the "Aquarium des Lagons" (Nouméa), which is 2 h distant from Bouraké, in individual hermetic zip bags (one for each colony of each species) containing seawater from the collection site. At the laboratory, fragments of A. tenuis and M. digitata were attached to nylon wires and suspended in two 200 L tanks, one for each site of collection (see below). Fragments of Porites sp. were mounted on labeled 2 × 2 cm PVC plates using epoxy resin (Holdfast, Aquarium Systems) and placed at the bottom of the tanks. The exposed skeleton was covered with the resin to avoid turf algae proliferation and potential skeletal dissolution. All 360 fragments (2 sites, 3 species, 15 colonies per species, 4 samples per colony, resulting in a base replication of n = 15 per treatment) were allowed to recover for 3 weeks in the two flow-through aquaria (see below for more details) settled at the same temperature (26.0°C ± 0.5°C) recorded in situ at the time of collection. Seawater pH and carbonate chemistry replicated mean values of Bouraké and reference sites (ca. 7.7 and 8.1 pH NBS units, respectively). Temperature and pH were kept constant using heaters and bubbling pure CO 2 in each tank. All tanks were connected to an IKS logger system (IKS, AquaStar, accuracy ±0.05 pH unit and 0.1°C) continuously monitoring but not recording the pH conditions.
During the recovery period, coral fragments received the same lighting intensity (ca. 250 μmol photons m 2 s −1 ) using four Aquablue Plus neon bulbs (15.000 K, Giesemann), and they were fed once a week with freshly hatched Artemia salina nauplii (regular equivalent concentrations tossed in the tanks; Houlbrèque et al., 2015). At the end of the 3 weeks recovery period, the naked skeleton of branching corals was already covered by the tissue, while new tissue was visible on the resin embedding Porites sp.

| Aquarium conditions and experimental setup
At the end of the recovery period, corals from each site of collection (called origin hereafter) were assigned to one of the 12 experimental tanks as follows. Three replicate 10 L tanks were set up for each of the four conditions of pH (monitored with an IKS logger system and pH Metrohm in NBS scale) ( Table 1) with averages recorded during the 100-day experiment: (i) stable pH (Control; pH NBS 8.11 ± 0.05; pCO 2 474 μatm); (ii) future reef pH based on the RCP7.0 IPCC scenario (IPCC Report, 2021) (Future; pH NBS 7.76 ± 0.07; pCO 2 1192 μatm); (iii) extreme stable pH (Extreme; pH NBS 7.54 ± 0.08; pCO 2 2115 μatm); and (iv) variable pH reproducing the semidiurnal tide phase variation of pH in Bouraké with four peaks of pH on a 24-h basis, two at low and two at high pH (Variable; mean pH NBS ranging from 7.56 ± 0.07 to 8.07 ± 0.07; pCO 2 ranging from 1968 to 533 μatm, respectively). While extreme stable pH was arbitrarily chosen, control and variable pH were settled based on the mean and variance data collected respectively at reference and Bouraké reef between 2016 and 2020 (see Maggioni et al., 2021). Five fragments per coral species (n = 3) and per origin (n = 2) were randomly positioned in each tank. Porites samples were positioned on the bottom of the tanks while the branching corals were suspended.
Seawater in the tanks was renewed at a rate of 16.5 L h −1 (renewal rate of ca. 165% h −1 ) and mixed using a submersible pump (micro-jet MC 320, Aquarium system). Each tank was supplied with seawater pumped from one to three 60 L sump tanks settled at different pH (7.4, 7.7, and 8.1). In addition, three more sump tanks were set up respectively at pH 7.6, 7.8, and 7.9 and were used to better simulate the pH variation in the variable treatment (see below). Each sump tank was continuously supplied with 50 μm-filtered seawater pumped at 5 m depth in front of the Aquarium des Lagons (Baie des Citrons, Nouméa). Each sump contained a submersible pump (23 W Eheim), a heater, and an external refrigerating system, both connected to an IKS logger system (IKS, AquaStar), and temperature and pH probes. Seawater pH NBS and temperature values were continuously monitored but not recorded by an IKS pH logger and daily verified (randomly during the day) using a portable pH-meter (Metrohm 826 coupled with a LL Aquatrode Plus SC, sensitivity of 0.001 pH units), both calibrated using NBS solutions (pH 4.0 and 7.0 from the Seawater National Bureau of Standards). Sumps were all maintained at 26.3°C ± 0.7°C, and at one of the 6 pH conditions, automated by IKS (AquaStar). Seawater pH was set to the desired value by bubbling pure CO 2 gas. Each experimental tank received seawater at the assigned temperature and pH from the respective sump tank. In contrast, for the variable pH treatment, the time function of the IKS system was used to simulate the pH changes measured at Bouraké, according to a semidiurnal tide cycle with high pH in the middle of the day/night (01:00-04:00 am/pm) and a low pH at the beginning of the day/night (08:00-11:00 am/pm). Briefly, tanks received seawater either from one or simultaneously from different sumps according to a predetermined timetable (Table S1; Figure S2). Irradiance in all tanks increased from 0 to ca. 250-300 μmol photons m −2 s −1 on a 12:12 h light:dark cycle (06:00-18:00 h lighting vs. 18:00-06:00 h darkness; ramp until 10:00 h and decline from 14:00 h) using LED lights (Mitras LX6100, GHL; see Table S2 for led brightness wave simulation) according to Biscéré et al. (2019), Houlbrèque et al. (2015 and Jacquemont et al. (2022). Irradiance was monitored during 48-h cycles using an NKE PAR with LI-193 spherical quantum sensor and showed maximum values ranging from ca. 270 to 320 μmol photons m −2 s −1 ( Figure S3). Tanks' positions were repeatedly changed during the duration of the experiment to minimize any potential spatial variability in the light intensity received by corals. A Seabird SeaFET pH logger and YSI 600 OMS-M probes were used to measure over 48-h cycles seawater temperature, pH T , salinity, and dissolved oxygen (DO) in the experimental tanks with the intention to document the variation in such a parameter over repeated diel cycles. Colonies were fed with freshly hatched Artemia salina nauplii once a week (regular equivalent concentrations tossed in the tanks) and were maintained under these experimental conditions for 100 days. Tanks, resins, and wires were regularly cleaned for turf algae proliferation.
Total alkalinity (A T ) was measured twice a month at various daytime to have the average total alkalinity reflecting the different pH values on each of the 12 tanks over the experiment. For that, water TA B L E 1 Seawater parameters measured during the 100-day experiment, and carbonate chemistry calculated for each pH condition (the three replicated tanks were pooled)  (2000), for KHSO 4 from Dickson (1990), for KHF from (Perez & Fraga, 1987), for [B] I from (Lee et al., 2010), and an averaged atmo- , NO x, and Si(OH) 4 were performed by colorimetry (Seal Analytical).

| Growth rate
After the recovery period and at the end of the experiment, each coral fragment was weighed using the buoyant weight technique (Davies, 1989). Samples were weighed using a Sartorius ENTRIS 224i-1S electronic balance (readability 0.1 mg) in seawater of known density (calculated from temperature and salinity, both measured in each 10 samples). Dry skeleton weight was calculated using the density of pure aragonite (2.94 g cm −3 ); growth rates (in terms of calcification rates) were calculated as the change in dry weight between the initial (T 0 ) and the final weight (T F ) during the experiment and expressed either in mg g −1 d −1 (normalized per weight of CaCO 3 ) or mg cm −2 d −1 (normalized per surface area of the coral, see below) depending on the species.

| Photosynthetic efficiency
During the last week of the 100-day experiment, maximum photochemical efficiency (F v /F m ) and the relative electron transport rate according to Ralph et al. (1999). Because corals were dark adapted we considered the first point the maximal quantum yield F v /F m . RLC uses a very short illumination period, which was found to be not correct to measure the real ETR (see Enríquez & Borowitzka, 2010).
Because we could not perform long RLC on the 135 corals, we limited our analysis to the maximum rETR max , assuming this limitation for all corals and conditions. Comparisons between treatments and origin were done on the maximum quantum yield F v /F m and the maximum rETR max .

| Photosynthesis and respiration rates
At the end of the experiment, seven corals for each pH condition (n = 4), each origin (n = 2), and each species (n = 3) were randomly selected (n = 168 total) and their oxygen production and consumption rates were measured daily in the light and dark during four consecutive days. Corals were individually placed in 100 mL Pyrex glass beakers which were filled with seawater from their respective tank treatment and hermetically sealed underwater with transparent cellophane and a rubber band (Jacquemont et al., 2022). Two control beakers without coral fragments were used to measure any metabolic activity of microbes in the water. We were able to process 28 samples simultaneously (plus two empty beakers) using two 15-place Rates of net photosynthesis (P net ) and respiration in the dark (R dark ) were calculated using the change in DO concentrations in each beaker corrected by the mean of the microbial activity measured in the two empty beakers, and normalized by the incubation duration (hours), the volume of seawater in each beaker (L), and the coral's surface (cm 2 ). Rates of gross photosynthesis (P g ) were calculated as: Data were normalized per surface area of the fragment as described below. P net , P g, and R dark are expressed in mg O 2 cm −2 h −1 and then converted in μmol O 2 cm −2 h −1 . Photosynthesis to respiration ratio (P g :R) was calculated using the value of daylight hours equal to 12 to calculate an integrated P g :R ratio for a 24 h day as follows: 2.2.5 | Tissue and surface measurements All fragments used to assess the coral photosynthesis and respiration rates were prepared and analyzed for their Symbiodiniaceae and chlorophyll contents. Protein measurements were performed on the same individual used during photosynthesis and respiration incubations for A. tenuis and M. digitata and from a new individual for Porites sp. Then, their skeleton's surface areas were measured. Coral tissue was extracted from the skeleton using an air pick in 20 mL filtered seawater and homogenized with a Potter tissue grinder.
For symbiont density measurement, 2 mL of the slurry was sampled to count the number of Symbiodiniaceae (n count = 8) using a Neubauer's cell under a stereomicroscope. 10 mL subsamples were centrifuged at 5000g for 10 min, the supernatant was discarded, and the pellet containing the symbiont was re-suspended in 10 mL of pure acetone to extract during 24 h at 4°C in darkness the chlorophyll a and c 2 . The solution was then centrifuged at 10,000g for 15 min and the supernatant was sampled to measure its absorbance at 630, 663, and 750 nm using a spectrophotometer (Evolution 201, Thermo Scientific). Chlorophyll a and c 2 concentrations were calculated using the spectrophotometric equations for dinoflagellates of Jeffrey and Humphrey (1975). Chlorophyll a and c 2 are given as total chlorophyll expressed in μg cm −2 .
Protein content was quantified using a BCA assay kit (Uptima, Interchim). Total protein was extracted according to (Hoogenboom et al., 2010)  The skeletal surface areas of samples were estimated using the paraffin wax-dipping method with two wax dips (Naumann et al., 2009;Stimson & Kinzie, 1991) for branching species (i.e., A. tenuis and M. digitata) and the aluminum foil technique (Marsh, 1967) for the massive species Porites sp.

| Symbiodiniaceae community analysis
2.3.1 | DNA extraction, PCR amplification, and sequencing Fragments of corals (n = 2-6 per species, depending on origin, and pH treatment) were collected and stored at −20°C. Total coral holobiont DNA (i.e., Symbiodiniaceae, polyp, and associated microorganisms DNAs) was extracted using a 2% CTAB-based protocol adapted from . The quantity and quality of extracted DNA were checked using a NanoDrop 2000 spectrophotometer (Thermo

| Statistical analyses and data presentation
Statistical analyses were conducted, and figures were produced using RStudio (R Development Core Team, version 4.1.0, 2021), including the packages ggplot2, ggpubr, car, and vegan. Homogeneity and the normality of variance distributions were tested using respectively the Levene test and the Shapiro-Wilkinson test and graphically verified with Q-Q plots. Statistical analyses were performed separately for each species as the three species we used are morphologically different. Each response variable was evaluated with separate linear mixed-effects models using the package lme4 (Bates et al., 2015). To evaluate differences in the main environmental parameters in the tanks during the experiment, i.e., pH, temperature, and nutrients measurements, the full model included pH as a fixed factor and tank as a random factor. To evaluate coral physiological responses, the full model included coral's origin (two levels: Bouraké and reference) and pH conditions (four levels: control, future, extreme, and variable) as fixed factors and tank as a random factor. The significance of fixed factors was evaluated with two-way ANOVAs (type III) using Satterthwaite's method. Significant differences between fixed factors were determined with Tukey's post hoc pairwise comparisons using the package emmeans. In contrast, F v /F m data did not meet the assumption of normality and they were compared

| Seawater parameters
Seawater pH, temperature, and carbonate chemistry were maintained at the target experimental values during the 100-day experiment (Table 1; Tables S3 and S4). There is a significant difference between conditions regarding temperatures (p < .001;

| Growth rate
Growth rates of the three coral species (Figure 2) were significantly different according to corals' origin (p ≤ .003) and pH conditions showing a p = .092 for pH (two-way ANOVA, Table 2 and Table S7 for all data analyses). The mean growth rates of the three corals originating from the reference and maintained at control pH, generally decreased when maintained at future, extreme, and also variable pH.
This decrease was more evident in the two branching corals than In contrast, the mean growth rates of the three corals originating from Bouraké and maintained at future, extreme, and variable pH (i) were always higher than their counterpart originating from the reference and maintained at the same pH condition (p ≤ .003); (ii) were almost always (with one exception, A. tenuis at variable pH) higher than their counterpart originating from the reference and maintained at the control pH condition; (iii) were almost always higher or did not change (with one exception, M. digitata at variable pH) when maintained at the control pH condition (Figure 2). For instance, the calcification rate of A. tenuis, Porites sp., M. digitata originating from Bouraké and maintained at the variable pH, which is what they normally experience at their site of origin, increased by 33.4% and 25.0% and decreased by 23.5%, respectively, in the control pH (pH NBS 8.11).

| Photosynthetic efficiency and electron transport rate
The F v /F m (i.e., maximum photochemical efficiency) and the rETR max (i.e., maximum relative electron transport rates) were significantly different for A. tenuis only (Table 2; Table S7). The F v /F m of A. tenuis differed according to both corals' origin (p = .001) and pH conditions (p = .018), whereas only corals' origin had an influence on rETR max (p < .001). Both rETR max (Figure 3) and F v /F m ( Figure S4) were higher for A. tenuis individuals from Bouraké, as compared to individuals F I G U R E 2 Reaction norm of the growth rates of corals from Bouraké (blue) and reference (gray) sites maintained during 100 days at four pH conditions (Control, pH NBS 8.11; Future, pH NBS 7.76; Extreme, pH NBS 7.54; and Variable, pH NBS 7.56-8.07). The arrow indicates the direction of the change from the condition of origin, assuming that corals from Bouraké might originate from seawater at either pH 7.7 (Future pH), 7.4 (Extreme pH), or fluctuating (Variable pH) and that corals from reference originate from seawater at pH 8.1 (Control pH). Data represented as dots are means of the three replicate tanks ± SE (overlapping SEs are in red) (n = 13-16, depending on species and pH condition; see Table S7 for all post hoc comparisons).
from the reference site, regardless of the pH treatment. The F v /F m for A. tenuis also showed higher values for future and extreme pH conditions than for control and variable pH conditions.

| Photosynthesis and respiration rates
Metabolic rates (i.e., P g , R dark , and P g :R) of the three coral species significantly varied between origins and pH conditions without a clear pattern of response. ANOVAs showed significant differences between origins in the P g (only Porites sp., p = .010), R dark (both A. tenuis, p = .022 and Porites sp., p = .009), but not for their ratio (p ≥ .355; Table 2; Figure S5; Table S7), and between pH conditions in the P g (p ≤ .008) and R dark (p ≤ .018) for A. tenuis and Porites sp. Overall, P g and R dark values were higher for future and extreme pH conditions for A. tenuis and Porites sp. (Figure 4; Table S7).

| Symbiodiniaceae, chlorophyll and protein content
A significant effect of pH (p = .044) was found only for Symbiodiniaceae density of the coral A. tenuis with higher density in the future condition and lower density in the variable condition (p < .050; Table 2; Table S7). Site of origin significantly affected the Symbiodiniaceae density and total chlorophyll content of A. tenuis (p ≤ .015; Table 2; Figure 5; Figure S6). Some interactions between pH and origin were also significant for A. tenuis (p = .035) regarding Symbiodiniaceae density and for Porites sp. regarding both Symbiodiniaceae density (p = .049) and total chlorophyll (p = .032).
Most of these differences were driven by the higher data variability in extreme condition, especially for corals originating from the reference site ( Figure 5; Figure S6; Table S7). Indeed, the post hoc com-
The pH treatments did not result in a significant change in the Symbiodiniaceae major type profiles for any coral species (p = .275; Table 3). However, an interesting observation was that in general corals from Bouraké had more consistent major ITS2-type profiles between replicate colonies, than corals from the reference site ( Figure 6). Individuals of the three coral species originating from the reference site had three times more type profiles when maintained at future or variable pH conditions than individuals under the control pH.

| DISCUSS ION
Corals from the Bouraké lagoon have evolved over generational time scales in an environmental condition (pH, temperature, and dissolved oxygen) chronically exceeding those predicted by the IPCC scenarios for the end of this century (IPCC Report, 2021).
The Bouraké geomorphology, coupled with an intense oxidative activity in the mangrove sediments, has been shown to cause seawater acidification, deoxygenation, and warming in the lagoon (Maggioni et al., 2021). Despite the extreme levels of acidification measured, which have been suggested to hinder coral calcification (Kleypas et al., 1999), an abundant, and well-diversified coral reef has developed in the lagoon. By assessing the holobiont physiological responses and the Symbiodiniaceae profiles of three coral species from both Bouraké and a reference reef to a large range in pH, we found that corals from Bouraké always exhibited higher growth rates and had a specific and more consistent ITS2 majority sequence than corals from the reference reef both under low and variable pH conditions. It seems likely that such patterns were linked to the strong life-long environmental fluctuations, which might have promoted coral resilience as previously suggested for other corals (e.g., Brown et al., 2022;Comeau et al., 2021;Enochs et al., 2020;Rivest et al., 2017;Schoepf et al., 2020). We are aware that this study explored only a few of the compensatory mechanisms that might be at the origin of the resilience observed for corals from Bouraké. We also admit that these coral populations may have developed such mechanisms in a much longer time scale than corals that will have to cope with rapid climate change. However, results from this study would suggest that adaptation to future OA conditions could be possible in the wild.

F I G U R E 3
Relative electron transport rate (rETR) measured during RLC curves of corals after 100 days of incubation at four pH conditions (Control, pH NBS 8.11; Future, pH NBS 7.76; Extreme, pH NBS 7.54; and Variable, pH NBS 7.56-8.07). Bouraké corals are in blue, while reference corals are in gray. Data are mean ± SE; n = 12-16, depending on species and pH condition; see Table S7 for all post hoc comparisons.

F I G U R E 4
Gross photosynthesis (P g ) and respiration (R dark ) rates of corals from Bouraké (in blue) and the reference (in gray) site at four pH conditions (Control,pH NBS 8.11;Future,pH NBS 7.76;Extreme,pH NBS 7.54;and Variable,). Data are median ± 25th and 75th percentiles (n = 7); see Table S7 for all post hoc comparisons.

F I G U R E 5
Symbiodiniaceae density and total protein content of corals from Bouraké (in blue) and the reference (in gray) site after 100 days of incubation at four pH conditions (Control,pH NBS 8.11;Future,pH NBS 7.76;Extreme,pH NBS 7.54;and Variable,).
Data are median ± 25th and 75th percentiles (n = 6-8, depending on species and pH condition). No data are available for the protein content of Porites sp. from Bouraké incubated at the Extreme pH condition. Asterisks indicate statistical significance at p < .05; see Table S7 for all post hoc comparisons.

| Coral potential resilience to ocean acidification: What do we learn from lifetime adapted corals?
In agreement with the consensus on the effect of OA on coral calcification , recently revised by Leung et al. (2022), we found that calcification rates of the both A. tenuis and Porites sp.
reference corals were significantly decreased when exposed to future and extreme pH levels. These results are consistent with coral species-specific responses to OA with respect to calcification. We observed species-specific responses to OA on corals from Bouraké although almost always (with one exception, A. tenuis at variable pH) showed higher calcification rates than their counterpart originating from the reference and maintained at the control pH condition. The marginal effect of OA observed in corals from Bouraké has been observed on some other coral species (e.g., Pocillopora damicornis; Porites cylindrica; Siderastrea siderea) in previous aquaria experiments (Aichelman et al., 2021;Bell et al., 2022;Brown et al., 2022;Comeau et al., 2013Comeau et al., , 2019 and has recently been the subject of meta-analysis and reviews (Bove et al., 2020;Leung et al., 2022).
Our data showed that all three coral species from Bouraké calcified from 19.0% (A. tenuis) to 67.3% (Porites sp.) more than their counterpart from the reference site, when maintained at control pH ( Figure 2a,b). This unexpected finding suggests that Bouraké corals have adapted to OA because they have been exposed to extreme conditions throughout their lives, and regain even higher rates of calcification than the same coral counterparts adapted to openwater pH conditions, once at actual open-water pH (i.e., ca. 8.10).
The innovative aspect of our study is that we compared corals adapted to ambient seawater pH with corals likely adapted, or at least fully acclimated, to an extreme, and fluctuating environment overlaying future climate scenarios. This study demonstrated that Bouraké corals have become more resilient to OA, likely using plastic and/or assimilated mechanisms such as changes in coral-associated microorganisms (Camp et al., 2020). Brown et al. (2022) incubated in aquaria with stable (218 ± 9 μatm) and variable daily pCO 2 amplitude (911 ± 31 μatm) conditions the coral Pocillopora damicornis from both a flat and a sloping reef, where mean daily pCO 2 amplitude differed (797 ± 20 and 399 ± 8 μatm d −1 , respectively). As in our study, they found higher rates of calcification for corals from the more variable environment. The authors measured lower secondary calcification (i.e., CaCO 3 density) on corals from the variable environment with less intracellular pH acidosis, as previously found for other corals (Comeau et al., 2021;Cornwall et al., 2018;Gibbin & Davy, 2014).
Clearly, this compensatory mechanism might have a cost. Among the potential hypotheses explaining how corals can cope with the additional energy required to maintain high calcification rates using for instance proton pumping (Guillermic et al., 2021;McCulloch et al., 2012), it is known that corals: (i) might boost their endosymbiotic algal production (Barott et al., 2015;Castillo et al., 2014;Schoepf et al., 2013), (ii) use their energy reserves, such as proteins and lipids Towle et al., 2015), and (iii) increase heterotrophy (Edmunds, 2011;Houlbrèque et al., 2015).
Our data do not consistently suggest that Bouraké corals have acquired a particular mechanism that accounts for better physiological plasticity to cope with low pH conditions, thereby maintaining higher calcification rates. Indeed, we did not find a clear effect of coral origin with a systematic increase in photosynthetic rates, higher contents of Symbiodiniaceae, chlorophylls, and protein in Bouraké corals. However, our study confirms previous findings on the positive effect of high pCO 2 on photosynthetic rates, which might boost endosymbiotic photosynthates production for the coral proton pumping mechanism. It is well known that under acidified conditions at the cellular level, the photosynthetic activity of algal symbionts in coral increases until it exceeds the maximum level of CO 2 consumed by the algae (e.g., Barott et al., 2015;Gattuso et al., 1999;Gibbin & Davy, 2014). Although CO 2 -induced photosynthetic fertilization of symbionts may be less effective for M. digitata, A. tenuis, and Porites sp. showed higher photosynthesis under OA scenarios compared to the control condition. This was accompanied by higher ETR max , F v / F m values, and better exploitation of lower light intensities for A. tenuis. This species also had the highest growth rates, which highlights the symbionts involvement in coral growth. In our study, we suggest that higher productivity was supported by an increased concentration in the Symbiodiniaceae density and/or chlorophylls, but not consistently among species highlighting the variable effect of OA on productivity such as previously reported (Anthony et al., 2008;Wall et al., 2014). Moreover, the respiration rates were higher under acid-   (Agostini et al., 2021). We recognize that protein content can only partially describe the change in the coral energy reserves, since lipids and carbohydrates were not measured.
In our study corals were fed once a week with Artemia salina nauplii and such heterotrophic inputs might have helped them (Cohen & Holcomb, 2009;Drenkard et al., 2013;Edmunds, 2011;Houlbrèque et al., 2015). Indeed, it has been observed that the artificial diet we used, although limited compared to previous studies, helps corals maintain appropriate energy expenditure and calcify under acidified conditions (Houlbrèque et al., 2015;Houlbrèque & Ferrier-Pagès, 2009). This is one of the potential limitations of most existing experiments in aquaria. Indeed, it is difficult to imitate the natural contribution of zooplankton to the diet of corals in an aquarium.
However, all individuals were fed in the same manner and the nutrient levels (i.e., NO x , PO 4 3− , and Si(OH) 4 ) measured in the aquaria were quite similar between tanks and throughout the experiment.
Thus, the diet could perhaps explain some of the apparent resistance of corals to OA, but not the different responses with regard to calcification between Bouraké and the reference site.

| To what extent do pH fluctuations improve the physiological performances of corals?
One of our hypotheses to explain the success of Bouraké corals was the potential positive effect that diurnal pH fluctuations could have on coral metabolism, as previously found for corals in GBR mangrove lagoons (Camp et al., 2019), St Vincent and the Grenadines CO 2 vents (Enochs et al., 2020), and for others incubated in mesocosms (Brown et al., 2022;Dufault et al., 2012). Our results show that short-term exposure (i.e. acclimation of reference corals during our experiment) to fluctuating pH has a negative effect on coral calcification as all coral species from the reference decreased their growth when incubated at variable pH ( Figure 2c). In contrast, Bouraké corals, which are acclimatized and/or adapted to local conditions, when exposed to variable pH calcified from 7.6% (A. tenuis) to 116.2% (M. digitata) more than their counterpart from the reference site. Furthermore, Bouraké corals incubated at future and extreme pH maintained these higher growth rates and even increased them when grown under control conditions, confirming their resistance to OA. Although the duration of our experiment was longer than most OA experiments (i.e., 77% lasted 1-11 weeks; Brown et al., 2022;Ziegler et al., 2021), we acknowledge that the duration of variable pH exposure experienced by corals from the reference site during the course of the experiment was too short to compare with that experienced by corals at Bouraké. Corals exposed to a variable environment throughout their life are physiologically more plastic than corals adapted to stable environments (Kenkel & Matz, 2017); such plasticity is probably time-dependent. Future experiments should take into account the length of exposure to variable pH (as well as other environmental parameters) that the corals experienced before the collection.

| Does their Symbiodiniaceae community raise the physiology of Bouraké corals?
It is thought that species-specific metabolic responses to environmental stress could be due to different symbiont communities hosted by corals (Barott et al., 2015;Ziegler et al., 2015). Our data suggest that seawater pH level, whether stable or variable, at future or extreme levels, does not affect the Symbiodiniaceae major ITS2type profiles of corals (at least over 100 days), since at the end of the experiment we found no significant effect among treatments.
Interestingly, however, Bouraké corals exhibited more consistent major ITS2-type profiles between replicate colonies at low and variable pH than corals from the reference site. Microbiome stability (both for symbionts and/or bacteria) was linked to greater physiological resilience to OA (Ge et al., 2021;Grottoli et al., 2018;Quigley et al., 2017Quigley et al., , 2019Ros et al., 2021), suggesting that consistent major ITS2-type profiles for Bouraké corals under pH treatments could facilitate their success; however, further work will be needed to verify this hypothesis.
We observed distinct major ITS2-type profiles between coral species and native habitat (e.g. Bouraké versus reference site).
All coral species in this study were associated with the genus  et al., 2003), and exhibit apparent capabilities to adapt to locally stressful and/or fluctuating environments (Howells et al., 2012;Ng & Ang, 2016;Schoepf et al., 2015). Corals from Bouraké exhibited lower Symbiodiniaceae density than the reference corals but had higher photosynthetic efficiency. Corals associated with Cladocopium goreaui (C1) have previously shown elevated photosynthetic efficiency with regard to other Symbiodiniaceae species (Cantin et al., 2009;Morgans et al., 2020;Stat, Morris, & Gates, 2008;Wall et al., 2020). It is therefore possible that elevated photosynthesis is a common trait within some taxa from the C1 radiation and may explain the higher photosynthetic efficiency mechanism of C1 to face OA as observed by Ge et al. (2021) in Acropora valida. The elevated photosynthetic capacity of A. tenuis in Bouraké, coupled with the dominance of major ITS2 sequence C1, could explain their significantly higher growth rates, as previously observed in A. tenuis juveniles (Cantin et al., 2009;Little et al., 2004). The C1 ITS2 sequence was still dominant for some colonies of A. tenuis from the reference site (albeit different ITS2 type profiles that suggest different species), but in most instances, ITS2 sequences from the C3 radiation (i.e., C3k and C3bo) were most abundant; these are not among the most efficient in photosynthetic capabilities (Hoadley et al., 2016). We note the uncommon prevalence of C15 in A. tenuis ( Figure S7). We are confident from our negative control gels that this is real ( Figure S8) (Camp et al., 2020) and from a mangrove lagoon on the Great Barrier Reef (Camp et al., 2019).
This is in agreement with the recent finding of Hoadley et al. (2021) demonstrating that within C15, diverged lineages exist for Porites sp. living in different reef habitats and that it brings physiological differences.
(iii) M. digitata from Bouraké had nine major ITS2 type profiles, with C15/C15vi-C15vj-C15f-C15he the most common. Cladocopium of the C15 radiation are considered physiologically resistant (Fisher et al., 2012;Fitt et al., 2009;Nitschke et al., 2018). The same C15 major type profile was also found for the reference corals, alongside C73-C73a-C21. Coral association with C73 is rarely found and appears to be present according to relative light levels (Stat, Loh, et al., 2008). With this in mind and according to the equivalent light irradiance measured between tanks in our experiments, we could assume that light regimes of the native environment of corals modified their Symbiodiniaceae community associations. For instance, turbid environments such as Bouraké have been observed to limit the diversity of the associated symbiont community and could lead to predominant local adaptive genotypes (Smith et al., 2020). While it is beyond the intent of this study to determine which environmental drivers shape symbiotic communities, our data supports the growing body of evidence (Camp et al., 2020;Howells et al., 2016;Ziegler et al., 2017) that corals from distinct environments often have unique symbiotic partners that could be crucial to support their survival.
While we cannot define the role of each Symbiodiniaceae ITS2profile in the coral stress response, the difference in symbiotic In conclusion, we show that under the OA scenarios tested here, corals adapted to an extreme and variable environment have systematically higher calcification rates than reef corals living at a more stable environment. Although not all physiological measurements revealed significant differences between coral's origin, Bouraké corals showed that their calcification rates have been somewhat boosted, regardless of pH level and variability. Concomitantly, they maintain divergent symbionts with more consistent Symbiodiniaceae communities than corals from a more stable environment. The comparison between corals that have been exposed during their life to variable pH with corals exposed only during the 100-day incubation suggests that this is not rapid mechanism. Upcoming studies should explore the effects of the exposure duration (i.e., short vs. long-term) to variable pH on coral biomineralization mechanisms, which already appear to differ between more stable and fluctuating environmental conditions of pH (Comeau et al., 2021). The different dynamics of symbiotic partners might greatly influence calcification processes via distinct translocation of photosynthates ultimately influencing mechanisms of internal pH regulation (Allen-Waller & Barott, 2023;Cameron et al., 2022;Venn et al., 2022). The use of isotopic tools such as the boron pH proxy in the coral skeleton would help that way, via tracing and defining the inputs of specific symbiotic partners into biomineralization mechanisms.

ACK N OWLED G M ENTS
We would like to thank Richard Farman, Sylvain Govan, Jeff Dubosc, and the technical staff of the Aquarium des Lagons (Nouméa) for their welcome in their facilities and for their assistance during our aquaria experiment. We wish to express our thanks to captains of the IRD vessels and the diving technical staff for their fieldwork assistance. We are indebted to Clarisse Majorel for her help during coral DNA extraction and amplifications for metagenomic analyses. We are grateful to Rafael Valente and Christine Sidobre for their help during the aquaria experiment.
We thank the Laboratoire des Moyens Analytiques LAMA (IRD, Nouméa) for nutrient analysis. We also thank the Province Sud of New Caledonia for sample collection permits (no. 3413-2019) and Greg and Esmé (Chez Esmé) for their support and hospitality during field missions.

FU N D I N G I N FO R M ATI O N
CT received a PhD fellowship from the University of Western

DATA AVA I L A B I L I T Y S TAT E M E N T
Additional supporting information may be found online in the Appendix S1 for this article. All data are available from the Dryad